Self-Consistent Approximations for Superconductivity 
beyond the Bardeen-Cooper-Schrieffer Theory 
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We develop a concise self-consistent perturbation expansion for superconductivity where all the 
pair processes are naturally incorporated without drawing "anomalous" Feynman diagrams. This 
simplification results from introducing an interaction vertex that is symmetric in the particle-hole 
indices besides the ordinary space-spin coordinates. The formalism automatically satisfies con- 
servation laws, includes the Luttinger-Ward theory as the normal-state limit, and reproduces the 
Bardeen-Cooper-Schrieffer theory as the lowest-order approximation. It enables us to study the 
thermodynamic, single-particle, two-particle, and dynamical properties of superconductors with 
competing fluctuations based on a single functional $[G] of Green's function G in the Nambu space. 
Specifically, we derive closed equations in the FLEX-S approximation, i.e., the fluctuation exchange 
approximation for superconductivity with all the pair processes, which contains extra terms besides 
those in the standard FLEX approximation. 
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I. INTRODUCTION 

The remarkable success of the mean-field Bardeen- 
Cooper-Schrieffer (BCS) theory in describing classic 
superconductors 1 ^ may be attributed to the fact that 
an infinitesimal attraction between electrons suffices to 
form Cooper pairs responsible for superconductivity. 3 
On the other hand, we now have an increasing num- 
ber of superconductors and Fermi supcrfluids where 
relevant interactions are stronger beyond the applica- 
bility of the mean-field theory. They include super- 
fluid 3 He at high pressure;^ heavy-fermion supercon- 
ductors with competing fluctuations^— high-T c cuprate 
superconductors^ and the BCS-BEC (Bose-Einstein con- 
densation) crossover in trapped atomic gasesJ^ The 
present paper is devoted to developing a convenient 
systematic method to describe those systems based on 
many-body quantum field theory. 

To be specific, we will focus on a two-body interac- 
tion and extend the Luttinger-Ward self-consistent per- 
turbation expansion for normal states^ to superconduc- 
tors in such a way that, given a normal-state approxi- 
mation, all the pair (i.e., "anomalous") processes deriv- 
able from it are automatically incorporated. Reproduc- 
ing the Hartree-Fock theory as the lowest-order approx- 
imation, the Luttinger-Ward formalism enables us to in- 
clude correlation effects systematically and microscopi- 
cally up to a desired order. Its key ingredient is a func- 
tional <£> [G] of the single-particle Green's function G de- 
scribed by closed skeleton Feynman diagrams, as prop- 
erly identified and called "^-derivable approximations" 
by Baym. 1 - Indeed, choosing an approximate $ enables 
us to calculate the whole thermodynamic series rang- 
ing from the thermodynamic potential and single-particle 
Green's function to the two-particle and higher-order 
Green's functions^ 3 - To put it another way, $ uniquely 
determines how to resolve the Bogoliubov-Born-Green- 
Kirkwood-Yvon (BBGKY) hierarchy- 3 '- 4 An additional 
prominent advantage of the formalism is that it automat- 



ically satisfies conservation laws, as shown by Bayrm 12 
Hence, one can also study non-equilibrium phenomena 
by only changing the imaginary-time Matsubara contour 
into the real-time Schwinger-Keldysh contour . 13 i 15 ~ — 

It should be noted that a formal extension of the 
Luttinger-Ward theory to superconductors was per- 
formed by de Dominicis and Martin already in 1964. 19 
However, they did not present any practical method of 
how to efficiently collect pair processes with anomalous 
Green's functions, whose number grows exponentially 
in comparison with normal diagrams as we proceed to 
higher orders. For example, (a)-(j) in Fig. Q] exhaust the 
third-order diagrams for $ in terms of the symmetrized 
vertex r^°^ of Abrikosov et al.^ where anomalous dia- 
grams (c)-(j) are already four times as large in number 
as normal diagrams (a) and (b). 

Here, we develop a perturbation expansion of gath- 
ering all the pair processes systematically and concisely 
without drawing anomalous diagrams for both the singlet 
and triplet pairings or their mixtures. This is made pos- 
sible by introducing a bare interaction vertex V} ^ that 
is further symmetrized in the particle- hole (i.e., creation- 
annihilation or "charge") indices. Thus, processes (a)-(j) 
in Fig. [1] can be represented by a single diagram on the 
right-hand side with a definite analytic expression, i.e., 




FIG. 1: Third-order diagrams for <f>. A square (a small circle) 
represents the bare vertex r' ' (T^ ') that is symmetric in the 
space-spin coordinates (both the space-spin coordinates and 
particle-hole indices). The number below each of diagrams 
(a)-(j) denotes its relative weight. 
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eq. (|26a[) below, which reproduces their weights correctly. 
The key functional $ is defined as a closed skeleton ex- 
pansion in terms of the matrix Green's function G in 
the Nambu space, which may be approximated by some 
selected terms from the whole series in practical calcu- 
lations. Besides the Dyson-Gor'kov equation for G and 
an expression of the thermodynamic potential, we also 
derive the Bethe-Salpeter equation for the two-particle 
Green's function that can be solved based solely on $. 
The formalism enables us to study not only supercon- 
ductivity itself but also competing effects of charge, spin, 
and superconducting fluctuations or orders simultane- 
ously through the matrix structure of T^°\ This last 
point may be regarded as a definite advantage of the 
present approach over the one based on the Stratonovich- 
Hubbard transformation, 21 where only a single fluctua- 
tion relevant to the auxiliary field may be taken into ac- 
count after some standard approximation. 

Using the formalism, we also extend the fluctuation 
exchange (FLEX) approximation originally developed 
for normal states 2 ^ - — so as to incorporate all the pair 
processes for the spin-singlet pairing. The FLEX ap- 
proximation has been successful in describing anoma- 
lous normal-state transport phenomena of high-T c super- 
conductors as well as some organic and heavy-fermion 
superconductors^ Its extension to superconductivity as 
described above, which has not been performed yet and 
will be called "FLEX-S approximation" below, will be 
useful for understanding those superconductors quanti- 
tatively. This will be carried out in a transparent and 
concise manner in terms of r*- -* in contrast to the origi- 
nal derivation for normal states 22r21 

Naturally, considerable efforts have been made to go 
beyond the BCS theory. Eliashberg incorporated the 
electron-phonon interaction explicitly as the source of 
the attraction to derive "strong-coupling" equations of 
superconductivity within the second-order perturbation 
in terms of the renormalized electron Green's function 
and bare phonon propagator. 2 - Leggett seriously consid- 
ered interactions between particles other than the pair- 
ing part to develop a theory of superfluid Fermi liquids 
at low temperatures^ Anderson and Brinkmart 2 ^ pre- 
sented the idea of "feedback effects," i.e., a change of the 
pairing interaction through the superfluid transition, to 
understand the A phase of superfluid 3 He that is realized 
against the mean-field Balian-Werthamer theory^ The 
idea was elaborated by Brinkman et aL 30 and Kuroda^i 
who both incorporated the first few anomalous processes 
into the normal-state pairing interaction composed of 
the particle-hole "paramagnon" diagrams™ and evalu- 
ated the extra free energy as a power series in the p- 
wave energy gaps. Subsequently, Tewordt^ 3 . refined this 
approach into a self-consistent ^-derivable approxima- 
tion, where he included those pair processes that can 
be obtained from the normal particle-hole diagrams by 
successively replacing a pair of normal Green's func- 
tions connecting adjacent vertices by a pair of anoma- 
lous Green's functions; relevant diagrams in the third 



order are (a), (c), (f), and (j) of Fig. Q] Tewordt 
and coworker a 34 ' 35 later applied the formalism to the 
two-dimensional Hubbard model to clarify quasiparticle 
and spin excitations of high-T c superconductors. The 
same formalism was used by other groups for the two- 
dimensional Hubbard model 3 ^— and d-p modeler— On 
the other hand, Eagles^ 2 - and Leggett^ 3 - discussed the 
BCS-BEC crossover of the isotropic spin-singlet pairing 
at zero temperature by combining the BCS gap equation 
with the equation for the chemical potential. The tran- 
sition temperature of the BCS-BEC crossover problem 
was calculated subsequently by Nozieres and Schmitt- 
Rink 4 ^ based on the Thouless criterion. 45 Haussmann 
and coworker a 46 ' 47 studied this problem with an im- 
proved ^-derivable approximation to include interaction 
effects among unpaired and paired fermions. Specifi- 
cally, their $ consists of the normal particle-particle lad- 
der diagrams appropriate for low-density systems, plus 
those anomalous processes obtained from the formers 
by successively replacing two pairs of normal Green's 
functions connecting adjacent vertices by two pairs of 
anomalous Green's functions; relevant diagrams in the 
third order are (b) and (h) of Fig. 1, disregarding (e) 
and (g) which may be important even within the ladder 
approximation. We should also mention an alternative 
approach on the issue where some effective pairing inter- 
action is used to solve the Eliashberg equations. 48 How- 
ever, there is an ambiguity as to how to construct the 
effective interaction. Similarly, the FLEX approximation 
has sometimes been augmented rather phenomenologi- 
cally to explain the pseudogap phenomena observed in 
underdoped cuprate superconductors ! 49 ' 50 To be specific, 
the "FLEX+T-matrix" approximation 5 ^ consists of (i) 
identifying the pairing interaction obtained by the FLEX 
approximation as the "bare" interaction in the Thouless 
criterion and (ii) incorporating the contribution of the 
two-particle "superconducting fluctuations" additionally 
into the single-particle self-energy. However, the proce- 
dure may contain some double counting of elementary 
processes and has yet to be examined on a firm micro- 
scopic basis. Note in this context that the Thouless crite- 
rion was originally derived by the particle-particle ladder 
approximation in the bare perturbation expansion for a 
two-particle Green's function, 45 reproducing the same T c 
equation as the mean-field BCS theory. Thus, one may 
be convinced that we still do not have practical meth- 
ods of incorporating all the anomalous processes that are 
naturally present below T c . 

This paper is organized as follows. In §2, we develop 
a basic formalism to calculate the thermodynamic po- 
tential, single-particle Green's function, and two-particle 
Green's function of superconductivity self-consistently 
based on a single functional $ [G] , whose expression may 
be obtained by a concise perturbation expansion in terms 
We carry it out in the coordinate representation 
so that general inhomogeneous systems may be handled. 
In §3, we focus on homogeneous systems and transform 
the results of §2 into the momentum- "energy" representa- 
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tion. Closed equations in the FLEX-S approximation are 
also derived. Section 4 presents a brief summary to indi- 
cate main results. The whole contents of the present pa- 
per may be regarded as the Fermi-superfluid counterpart 
of the self-consistent perturbation expansion recently de- 
veloped for BEC— ~— We set % = ks = 1 throughout 
with fee the Boltzmann constant. 



II. BASIC FORMALISM 

The system we consider consists of identical particles 
with mass m and spin 1/2 described by the Hamiltonian 



with 



H — H + H int , 



H = J dtitfi&Knl/fa), 



(1) 



(2a) 



Hi, 



d£i / d^V' t (ei)^ t (6)^(r 1 -r 2 )^(6)V(ei). 

(2b) 

Here £i = (ri,a\) with r% and oti denoting the space 
and spin coordinates, respectively^ tp* and tp are the 
creation and annihilation operators of the fermion field, 
respectively, — K\ = —Ti 2 V\/2m — [i with \x the chemical 
potential, and V is the interaction potential. Though dis- 
regarded here, the effect of a lattice or trap potential may 
be included easily in K\ . The formulation in this section 
will be carried out in the coordinate representation so 
that general inhomogeneous systems can be handled. It 
is also applicable to lattice models such as the Hubbard 
model, for which every integration over r above should 
be replaced by a summation over the lattice sites. 

Let us introduce the Heisenberg representations of the 
field operators in the Matsubara formalism by22, 



(1) = e^ H V(a)e- 



^2 H (1) 



^(£l)< 



■TlH 
-Tiff 



(3) 



where argument 1 in the round brackets is defined by 
1 = (£l,ti) and variable ti lies inO<ri</3=l/T with 
T denoting the temperature. The operators an d 
?pf(l) are denoted as ip(l) and ^(1) by Abrikosov et al. 20 , 
respectively; distinguishing them by the subscript i = 1, 2 
(i.e., the particle- hole or "charge" index) as above enables 
us to simplify the notation and formulation considerably, 
as seen below. 



A. Green's function, Dyson-Gor'kov equation, and 
thermodynamic potential 

Using ^P(l), we define Green's functions by 

G y (l,2) = -(T T ^ H (l)V3 H -i(2)>, (4) 



where T T is the "time" -ordering operator and (■ • • ) de- 
notes the grand-canonical average in terms of The 
standard normal-state Green's function G corresponds to 
Gii. The elements satisfy 

Gy(l,2) = G 3 ,, ,;2. 1: = G^fon, £it 2 ), (5) 

with the superscript * denoting the complex conjugate. 
The Nambu matrix G = (Gy ) is written explicitly as 



G(l,2) 



G u (l,2) G 12 (l,2)' 
G 2 i(l,2) G 22 (l,2) 



(0) 



It obeys the Dyson-Gor'kov equation 

[Gq X (l, 3) - £(1, 3)] G(3, 2) = <x 5(l, 2), (7) 

where <To is the 2x2 unit matrix, 5(1, 2) = 5(ti — t 2 )5(£i — 
£2)) and integrations over barred arguments are implied. 
The quantity Gq 1 is defined by 



<V(1,2) 



where (73 is the third Pauli matrix. Finally, £ in eq. ([7]) 
denotes the self-energy matrix 



£(1,2) 



E11 (1,2) E 12 (l,2)' 
E 21 (l,2) £22(1,2) 



(9) 



It follows from eqs. ©, ©, and (|5J) that £ satisfies the 
same symmetry relations as G in eq. ([5]). 

As shown by de Dominicis and Martin^ we can ex- 
press the thermodynamic potential ft = — fi^ 1 lnTre~^ ff 
as a functional of G, i.e., f2 = £l[G], so as to satisfy 



sn 



= 0. 



(10) 



Let us write Q[G] as 



O = -^Tr[ln(-G - 1 + E)+EG] (11) 

where £ = £[G], $ = $[G], and the operator Tr for the 
Nambu matrices is defined by 



TTA = A n (l,l+) + A 22 (l + ,l) 



(12) 



with the subscript of 1+ denoting an extra infinitesimal 
positive constant in t\ to place creation operators to the 
left of annihilation ones for the equal-time average. Using 
eqs. (Jjj and (fTT|) . one may show easily that eq. (fTO)) is 
transformed into a relation between E and $ as 



£„(l,2) = 2/3 



5$ 



^(2,1) 



(13) 



Equation (|lll) with eqs. ([7]) and (IT31 forms an exten- 
sion of the normal-state Luttinger-Ward functional^ to 
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FIG. 2: Bare interaction vertex r' ' that is symmetric in a 
pair of outgoing (incoming) lines. 



superconductors. Note that we have given as a func- 
tional of G following Baym 12 instead of E in the original 
treatment 

It follows from eqs. (J7J, (|TT|) . and (fT5)l that we can 
calculate G and f2 self-consistently once $ = $[G] is 
given explicitly. The particle number N is obtained by 



N = ^Tr Ga 3 , 



(14) 



which may be used to change an independent variable 
from fi to N by the Legendre transformation F = Q+/a/V. 



B. S-matrix and functional $ 

Perturbation expansions of SI and G in terms of Hint 
in eq. (|2b[) can be carried out conveniently with the S- 
matri3*2£ 



S = T T exp 



-^^(ll',^')^®^®^^)^^) 



(15) 

Here ^i(l) (i = 1,2) are the interaction representa- 
tions of the field operators given explicitly by ipi(l) = 
e nHo ^(£i)e~ nffo and ^ 2 (1) = e Tlff °V> t (£i) e ~ Tli?0 ' and 
F^ denotes the symmetrized bare vertes^S 

r^(ll',22') = ^(r 1 -r 2 ) ( 5(T 1 -r 2 ) 

x[5(l,l')5(2,2')-5(l,2')5(2,l')], (16) 

satisfying r<°) (11', 22') = rW (22', 11') = rC°) (l'l, 2'2) = 
-r(°)(12',21'). The interaction in eq. (HSJ can be ex- 
pressed graphically as Fig. [5J 

Now, the key functional <£> in eq. (fPTj) is given exactly 
in terms of eq. (|T5"|) as 



$[G] 



< S > 0c -1 







(17) 



skeleton, Go — ^C? 



where subscript o denotes the thermodynamic average 
with Hq and another subscript c implies retaining only 
connected Feynman diagrams. The right-hand side may 
be expressed graphically by the skeleton diagrams (i.e., 
diagrams without self-energy insertions) in the perturba- 
tion expansion for Q, with Go replaced by G.— Approx- 
imating eq. ()17[) by a few terms or some partial series, 
we can construct a self-consistent approximation, as al- 
ready noted in the last paragraph of mi Al Equation (fl7|) 
was proved by Luttinger and Ward for normal states . 11 ' 13 
However, it also holds true for superconductors by includ- 
ing all the pair processes. 



CO 



(a) 2 



(b)l 



FIG. 3: First-order diagrams for $. The number below each 
diagram denotes its relative weight. 



The Feynman rules for the perturbation expansion of 
eq. ([TT]) with eq. (| 15[) was given by Abrikosov et al^ 
for normal states. Following them, however, we need to 
calculate numerical factors from a separate expansion in 
terms of V. This may be the reason why the expansion 
with r(°) has not been used widely. It is shown in Ap- 
pendbJAl that the rules can be simplified considerably so 
that one may perform it directly in terms of T^-°\ even 
for superconductors. For example, the first-order contri- 
bution to $ is graphically given by Fig. [31 which can be 
expressed analytically as 

*i = ^rW(ll',22')[2G(l',l)G(2',2) 

+^(T,2)F(1',2')], (18) 

where G, F, and F are defined together with G by 
_ G n (l,2)-G 22 (2,l) 



G(l,2) = 



_ G 12 (l,2)-G 12 (2,l) 
F(l,2) = , 

- M ^ G 21 (l,2)-G 21 (2,l) 



(3(2,1), (19a) 
(19b) 
(19c) 



With eq. (|19[) . we have incorporated the first symmetry 
of eq. (|5|) manifestly into $i for the purpose of deriving 
r (0) below. 



C. Concise perturbation expansion for $ 

The expansion with is still cumbersome for collect- 
ing anomalous diagrams in superconductivity, however, 
as seen from Fig. [1] Hence, we further modify the ex- 
pression of eq. (|15[) so that the asymmetry between the 
incoming and outgoing lines in Fig. [5] is removed. To be 
specific, we introduce a bare vertex T^) ■■, (11', 22') from 



eq. |T8]) by 



»*' ,jj 



(11', 22') = 2/3 



5G^(1',1)G^(2',2) 



(20a) 



A straightforward calculation using the symmetry of eq. 
(jTHJ) shows that results of all the above differentiations 
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FIG. 4: Bare interaction vertex T$ . ., (11', 22') that is sym- 
metric in four external lines. 
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(1) 





(4a) 



(4b) 



FIG. 5: First- to forth-order diagrams for $. 



may be summarized as 

r^ i ,(ii',220 = ^^ -i (; r-'ci i'. 22'. 



2 

<5 ; 



, 0jj ,\ 

5 U '5 jf r(°\ll',2'2) 



^%,^r(°)(12,l'2'), (20b) 



which satisfies 



r^,(ll',22') =rW^(22',ll') = -1^,(12', 21') 



= -^,3-^% 22') = -r^ J ., ii , )ii3 _ i (2'l' ) 21). (21) 
Now, we can express eq. (|15p alternatively as 



(0) 



S = T T exp 



1 



-— ri ? „,(ll', 22') 



xA^-,(l')V> 3 -«(l)^(2')V 3 -j(2) 



(22) 



where A/" denotes the normal-ordering operator of placing 
creation operators to the left of annihilation operators 
with a sign change per every permutation of a pair of 
adjacent field operators.— The equivalence between eqs. 
([15]) and (|22|) may be checked easily by substituting eq. 
(|20bj) into the latter and writing the resultant expression 
without the N operator. Note that M is only relevant to 
equal-time averages of the first order in the perturbation 
expansion. 

The interaction in eq. (|22l) may be expressed graphi- 
cally as Fig. [4j which is symmetric in the four external 
lines. Using eq. (|22j). we can perform the perturbation 
expansion for $ (and also for any other quantities) con- 
cisely in such a way that all the anomalous processes are 
incorporated automatically. Figure [5] enumerates first- 
to forth-order diagrams for Thus, each of the first- 
to third-order contributions is exhausted by a single dia- 
gram, showing manifestly the advantage of using eq. ([21} . 
Moreover, the symmetry of eq. (f2~Tj) brings the simplifica- 
tion that every possible Wick decomposition for a distinct 
diagram yields the same contribution. Hence, we only 
need to consider a convenient decomposition (i.e., con- 
nection of vertices) for each distinct diagram and mul- 
tiply the result by the number of possible connections, 
which may be calculated easily based on a combinatorial 
consideration. 

Now, the Feynman rules for the expansion of eq. (fl7|) 
with eq. (|22|) are summarized as follows. 



(a) Draw all the nth-order closed skeleton diagrams 
that are topologically distinct. For each such di- 
agram, associate the factor (— l) ,l+1 /n!12 n /3. 



(b) For each small circle, associate rl^-^AA', vv'). 

(c) Identify the number C na of possible connections of 
vertices for diagram a under consideration. 

(d) Consider a specific connection of vertices for dia- 
gram a where every primed argument is linked with 
an unprimed argument, and associate Gj,j(X', fj) for 
each line connecting A^, and fjj. 

(e) Identify the number £ na of permutations to realize 
the connection. With the choice in (d), we may find 
it by replacing every vertex as T^) ■- ■, (11', 22') — > 

6 iX 6 i >i6 j i8j,i6il, l')<5(2, 2')V(ri - r 2 )<5(n - n) and 
equating £ na with the number of closed particle 
loops in the resultant normal diagram with V ^ 

(f) Multiply the expression by (— X) lna C na . 

(g) When calculating the two-particle irreducible ver- 
tex r_( ir ) by eq. (IMl) below, replace every Green's 
function as 



G,j(l,2)->G y (l,2) = 



_ G ij (l,2)-G 3 _ i , 3 _ i (2,l) 



(23) 



so as to incorporate the first symmetry of eq. ([5]) 
manifestly in T^ 11 -* . 

Rule (g) is relevant only to two-particle and higher-order 
Green's functions. 

As an example of rule (c) above, let us consider dia- 
gram (3) in Fig. [3] Its Wick decomposition may proceed 
as follows: Start from an arbitrary operator ipi (labelled 
a\ operator in H?u) and find its partner (labelled b\ op- 
erator in HP Dt ) from 2x4 possible candidates; to connect 
another line between H^ nt and H^ nt , pick out a pair of op- 
erators a2 and 62 from 3x3 possible choices; to connect 
a pair of lines between H^ nt and H? nt , find the partners of 
63 and 64 (labelled c\ and C2 operators, respectively) from 
4x3 possible choices in H? nt ; finally, select one of the 2 
possibilities to connect a pair of lines between (03,04) and 
(03,04). Thus, the number C 3 is obtained as 

C 3 = (2 • 4 • 3 2 ) • (4 • 3) • 2 = 2!4 2 3 3 2. 

Now, rules (a)-(f) enable us to write down an analytic 
expression for <J> 3 as 



= 4 TT ^i\:::,,;ii^22';r; /i >2'.3i^,(3'.2 j 

xlf ) m (33', 44')G r , rn (4', 5)0^(5', 4) 



kk',U 

xF^L'nn' (55', 66')G fi , ? (6', l)G;, s (l', 6), (24) 
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with £3 = 3 from rule (e) for this connection. An ele- 
mentary calculation using eq. (|20b[) shows that eq. (pM)) 
correctly reproduces the weights of diagrams (a)-(j) in 
Fig. [TJ It is convenient at this stage to introduce matri- 
ces T (0) and x (0) by 



(ii^|r(°)|22' )=r£ .,,. (n',2'2) 



(25a) 



(ll^|x (0) |22; y > = -0^(1,2)0^(2',!'). (25b) 
Using them, we can express eq. ([24| concisely as 

(26a) 



$3 = ^Tr(r(V) 3 



where the basis of Tr is given by the bracket vectors in 
eq. (g5]| 

The above consideration may be extended to the other 
diagrams in Fig. [SJ where the numbers of rule (c) are 
identified as d = 3, G 2 = 4!, C Aa = 3!4 3 3 4 2, and C ib = 
( 4 G 2 ) 4 3 ■ 2 2 2 3 = 3 5 2 9 = 3!4 4 3 4 . We thereby obtain 



$i = ^Trr(°)x (0) , 



$ 2 = -— Tr (r(°) X (°)) 2 
12/3 y ~ - > 



$4 a = -^Tr(r(°) X (0) ) 



(26b) 



(26c) 



(26d) 



1 



xr^.«n'(55 / ,66 / )rffi. - s ,(7 , 7',88') 
y 



' pp m 

xG n (2', 3)G^(3', 2)G n , p (W, 7)G ra (7', 6) 
xG ?g -(l',8)G rr (8',4)G P „- l (4',5)G^(5',T). (26e) 



The equivalence between eqs. (|T8|) and (|26bp may be 
checked easily by using eqs. (fl"9j) . (|20bj) . and (f25|) . Now, 
one may be convinced that eq. (|22[) considerably simpli- 
fies the perturbation expansion of superconductivity. 



D. Self-consistent approximations 

As already mentioned, eq. (1171) forms a convenient ba- 
sis for systematic approximations for superconductivity. 
To be specific, Green's function and the self-energy are 
determined self-consistently by eqs. ([7]) and (TT31 based 
on some approximate <£>. The lowest-order approxima- 
tion corresponds to the choice $ « $1 from the series of 
eq. (|26|) . which is exactly the mean- field BCS theory.^ 
It may be improved by including the next-order term as 
$ s» $1 + $2, where two-body quasiparticle scatterings 
are included. Hence, adopting this approximation on the 
Schwinger-Kcldysh contour ) 13 ' 17 i 18 we can describe ther- 
malization of superconductors, for example. 



We now introduce the FLEX-S approximation, which 
adds the sort of diagrams (3) and (4a) in Fig. 5 up to 
infinite order besides $1 and $2- As seen from Fig. 1, 
it contains the normal particle-hole and particle-particle 
diagrams, i.e., diagrams (a) and (b), respectively, where 
every exchange process is automatically incorporated 
through r(°), plus all the pair processes derivable from 
them by successively changing directions of a pair of in- 
coming and outgoing arrows at a vertex, i.e., diagrams 
(c)-(j). Generalizing the consideration that led to eqs. 
(|26a| and (|26dj) . we find the nth-order contribution <fr na 
(n > 4) as 



( i>^ = ^45^^(r (0) x (0) r 



2n/3 



(27) 



The functional $ of our FLEX-S approximation is defined 
b y $ FLEX-S = $ i + $ 2 + $3 + E^=4 $ «a, i-e., 

1 ^ „cni rn\ 1 
IP 

"in(i + r(°)x (0) )-r (0) x (0) l ■ (28) 



*FLEX-S = ^TrF(V 0) + ^Tr(r(V 0) ) 2 



2/3 



The corresponding self-energy is obtained by eq. (|13[) 
with $ ss ^FLEX-S' wnere the differentiation may also 
be carried out graphically by removing a single line from 
the diagrams for <E> in all possible ways. Using eq. (|2"T]) . 
we arrive at the expression 



£ij(l,2) 



;i3 lS |r( )|24 7r )+-(i3^|r( ) x ( )r(°)|24 7r ) 



-2(i%|r(°)x^r( 1 )|24 J , 



Gk(3,4), 



(29) 



(30) 



where r'- 1 -' is defined by 

r« ee r (0) (i + x (0) r (0) ) _1 . 

E. Two-particle Green's function 



Next, we consider the two-particle Green's function 

K iV , jr (U',22') ee (T^F(l)^ H _ l ,(l')^ H (2)^(2')) 

-G^(M')G jy (2,2'), (31) 

whose poles define collective modes such as density, spin, 
and superconducting fluctuations. Functional $[(?] also 
enables us to calculate eq. (l3Tj) unambiguously. To see 
this, let us apply an artificial nonlocal potential W(l,2) 
described by the S-matrix 



S u ee T T exp 



-^(2)V^,(2')^,j(2',2) 



(32) 



Green's function in the presence of U = (Uij) is defined 
by22 



(T T S u j>f(l)^_ ¥ (V)) 
{S u ) 



(33) 



which reduces to eq. (HJ) as U — !• 6. Now, one may see 
easily that eq. (|3"Tj) is obtained from eq. ([55]) by 



/C M (ll',22') = -2 



5U rj (2',2) 



(34) 



This expression tells us that we only need to know the 
linear response of G u to W for obtaining the two-particle 
Green's function. 

To find 5G U , we start from the Dyson-Gor'kov equa- 
tion, which is modified from eq. ([7]) into 52 ' 54 ! 55 

[do Hh 2) - Z?(l, 2) - E"(l, 2)] G"(2, l') = a S(l, l'), 

(35) 



with 



£^,(1,1) = . (36) 



Let us change U —> U + 8U in eq. (|35|) and subsequently 
set U = 0. We thereby obtain the first-order equation 

G- 1 (l,2)8G u (2,l') = [5U' (1,2) + 5t u (l,2)]G(2,l'). 

(37) 

Using eq. (HUl) for the self-energy, we can express 5T, U 
above in terms of 5G U as 

«5E^, (1, 1') = (11', 22')<5G^(2', 2), (38) 



where r (11', 22') denotes the "irreducible" vertex 



-.(ir) 
ii' ,. 

defined by 



-.(ir)_ 



(11', 22') = 2(3 



5G iH (l',l)5G fj (2';2) 



(39) 



This differentiation should be performed after the re- 
placement of eq. (f2"3"|) in as already mentioned. Let 
us substitute eq. (f3"5]) into eq. (I3T1) . multiply the resul- 
tant equation by G(3, 1) from the left, and integrate over 
1. Changing the arguments appropriately, we obtain a 
closed equation for SG U as 



*G&(1, 10 = %(1,2)G^,(2', 10«/4j,(2,20 + %(1,2) 
x G^, (2', 1')^ ,,,(22', 3'3)<5G^, (3, 3'). 



(40) 



At this stage, it is convenient to introduce the nota- 
tions 



(ll^G" = 5G%(1,1'), (ll' tl ,\SU = 6Uu>(l,l% 



(41) 



and 



(llj i ,|£|22; j .,) = ^ i (ll , ) 2'2) ) (42a) 

(11^11^122^) ^r|^(ll',2'2), (42b) 
(ll^|x (0e) |22^)=-G«(l ) 2)G^(2MO 

+G 4 , 3 - J '(1,2')G 3 _^(2,1'), (42c) 

(11^1122^,)= 5^5(1,2)5(1', 2'). (42d) 




+ 




FIG. 6: Bethe-Salpeter equation for superconductivity. A 
black (shaded) circle denotes F (£/ ir '), and an internal line 
represents G. 



It follows from eqs. ©, ([SI]), dSSJ), and @2) that matrices 
M =jC, T (ir) , x (0c) all satisfy 



(ll^|M|22;. y ) = (2'2^|M|l'li' < ) 

= -(l'l 3 ^, 3 _ 4 |M|22' r; 



(43) 



The first equality is also obeyed by x^ of eq. (|25b[) . 

Using eqs. (|25bl) . ([56]). (|4"l~j) . and ^2j, we can express 
eq. gO} as (SG W = -±x (0o) ^ - x (0) r (ir) <5G w , or equiva- 
lently, 

5G U = -i(i + x (0) r (ir) rV 0e ^i7. 



Finally, we obtain an expression of /C by eq. (|34l) as 

(i + x (0) r (ir O"V 0o) . (44) 

It may be useful to introduce the full vertex 

r = r^(i + x (0) r (ir) )"\ (45) 

which is nothing but the Bethe-Salpeter equation^ for 
superconductivity and given graphically by Fig.Hl Using 
r, we can express eq. (|4"4l alternatively as 



!C = x {0c) -x (0) Lx {0c) - 



(46) 



Equation dm or 0BJ with eqs. (|25b| . ([3^1) . and @2) tells 
us that we can also calculate the two-particle Green's 
function of eq. ([3T]) once $ is given explicitly. 

It may be illuminating to write down r' lr ) explicitly 
in the FLEX-S approximation. Let us make the replace- 
ment of eq. ([23jl in eq. (|28|) . substitute the resultant $ 
into eq. ([39jl . and perform the differentiation. Also using 
the symmetry noted around eq. (|43[) . we obtain 



(ii^|r^|22^,) 
= (ii^|r(°)|22; y ) - 2[<i2 ij |r(°)xWr(°)|i% 

-<i^8-j'ir (0 ¥ 0) rW|i%,a-j>] 
+<i2 iJ -|r(°)x (0) r( 1 )|i'%,) 

-{l^ir^x^ll'^a-i) 
+2[<13 iS |rW|24'-,><2'4 j , r |r«|l'3^ / > 



-(13 l ,|r^^|2'4;_ j , r ,)(243_ jJ |r^^|l'3^, 
x(34' gr ,|x ( °)|3'4,, r ), 



(47) 



which is different from of eq. (|30|) . Thus, a clear 
distinction among r^, JZ , and T is necessary. 
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F. Alternative approach 

Before closing this section, we comment on an alterna- 
tive self-consistent approach adopted by Haussmann et 
al£L of expressing ft in terms of a renormalized vertex T 
besides G. Its theoretical basis was given by de Dominicis 
and Martin^ see also refs. l57l andl58l on this point. How- 
ever, the vertices T n obtained by de Dominicis^^ and 
used by Haussmann et alM- have no internal degrees of 
freedom describing only density fluctuations, in contrast 
toourr (0) of eq. (|3DJ|. Indeed, constructing ft = Q[G,T n ] 
with a matrix structure in r° seems to be a non-trivial 
problem. It will be practically more convenient to use 
ft = ft [G] that can easily be written down solely in terms 
of the smallest expansion unit G. Note in this context 
that: (i) the results obtained by using ft = ft[G, T ] 
should necessarily be reproduced by the approach based 
on ft = ft[G], as in the case of the particle-particle lad- 
der approximation by Haussmann et al.^ (ii) the vertex 
r a in ft is generally different from both the irreducible 
vertex T (ir) of eq. dSHJ) and the full vertex T of eq. (14"51 . 



III. HOMOGENEOUS SYSTEMS 

We now simplify the formalism of §2 for homogeneous 
systems. Lattice models may be handled similarly with 
modifications given in AppendbfBl 

A. Momentum- "energy" representation 

Equations (@| and © and the delta function 5(1,2) 
can be expanded for homogeneous systems as 

Gy(l, 2) = £ G ta ^(p) e<*C*-*>, (48a) 
p 

E tf (l, 2) = £ £ ia , i/3 (p) e^-^, (48b) 
p 

8(1,2) =8 a ^e^-^. (48c) 
p 

Here T\ - (rx, in) and p = (p, ie n ) with e n = (2n + 1)ttT 
the fermion Matsubara frequency (n = 0, ±1, • • • ), sub- 
scripts Q;( 3 denote spin indices, and the summation over 
p is defined by 



d 3 p 
(2^F' 



(49) 



We call the imaginary quantity ie n "energy" below. Sub- 
stituting eq. (l48l) into eq. 0, we can transform the 
Dyson-Gor'kov equation into the 4x4 matrix equation 

[Gq 1 (p)-£(p)]G(p) = 1, (50) 
where Gq 1 (p) and 1 are defined by 

(ia\do l @)\j0) = SijSafilien - (-l)'" 1 ^] (51) 



with e p = p 2 /2m—fi and (ia\l\j(3) = 8ij8 a p, respectively. 
The standard arrangement of the basis vectors for the 
direct product C x S of the particle-hole (i.e., charge) 
space C (i = 1, 2) and spin space S (a =\,\) is given by 
|csi) - |lt), |cs 2 ) = \U), |cs 3 > - |2t), and |cs 4 ) = |2|). 
Green's function in this representation reads 



G(p) = 



G(p) F(p) 

-Bp) -g(p). 



(52) 



where each quantity with an underline is a 2 x 2 matrix 
in spin space given as a Fourier coefficient of eq. (|19[) . 
It follows from eqs. © and (TT9"j) that these submatrices 
satisfy 



G a p(p) = G 0a (-p) = G} a (p*), 



(53a) 



F a p{p) = -F Pa {-p) = -FZJF). (53b) 
The self-energy matrix may also be expressed as 



Up) Mp) 
-Hp) -Up). 



(54) 



whose submatrices clearly obey the relations of eq. (|53|) . 

Next, we transform £ = £[G] into the momentum- 
energy representation. Since the expression depends on 
the approximation we adopt, we specifically consider the 
FLEX-S approximation given by eq. (|29|) . To begin with, 
we expand V(ri — r2)8(n — T2) in eq. (fT6|) as 



V(ri - r 2 )6(n - r a ) = £ V q e***" 



(55) 



where q = (q, iui) with uji = 2£wT the boson Matsub- 
ara frequency (I = 0,±1, • • • ), and the summation over 
q is defined in the same way as eq. ([4TJ]) with the re- 
placement p — » q and n — > t. Let us substitute eq. (jTHJ) 
into eq. (|20b[) . use eqs. (|48c[) and (|55| subsequently, and 
make some changes of integration variables; see also Ap- 
pendlxlBl for details. We thereby obtain an expansion of 
eq. (|25a|) as 



(ll^|r (0) |22^.,) = "£{ia,i'a'\T( \p,p',q)\jPJ'f3') 



pp'q 

Y p if-fx — if- — if -r 2 +if -Pa 



where p- is defined by 

P- = (p- q,ie n - iwt), 
and the Fourier coefficient is given by 



(56) 
(57) 



(ia,i'a'\T^{p,p',q)mf0') 
[(-l) i+i 8 H ,6 jj ,5 a , a 8f l ,f i V q 



1 



— (— SijSi' ji S a j3S a ' p'Vp-pi 



+(— 6ji£-i8j,a-ii6 a pi6 a tpVp+p>-q\ 



(58) 



9 



Note that the arguments t\, t[, t 2 , and t' 2 in eq. (|56)l 
are properly expanded in fermion Matsubara frequencies; 
this is why we have expressed eq. (|55|) in terms of bo- 
son Matsubara frequencies. Using eq. (|48a[) , we can also 
transform eq. (|25b[) into 



pp'q 



xe 



ip-r± — -r 1 — ip -r-2+ip _-r 2 



(59) 



with 



(ia,i'a'\x (0) (P,<l)\3(3,fP') = -G ia ,jpip )G rp ^, a , (p_). 

(60) 

Substituting eqs. ([15]). and dSSJ) into eq. P9"]h we 

obtain the self-energy for homogeneous systems in the 
FLEX-S approximation as 



r (0) (p )P ,q) 



+ 1 E E (0) (p> pi > 9k (0) (ft . ^E (0) (pi . p> D 
pi 

-2 £ r(°) (p, Pl , 9 )x (0) (pi , ( Pl , p, q) 

Pl 



(61) 



where F^ is the momentum-energy representation of eq. 
([50)1 satisfying 



E (1) (p, p',g) = r(°> (p, p',g)-J2 £ (0) fa pi > 9) 



pi 



xx (0) (pi,9)r( 1 )(p 1 ,p', q ), 

and c*x (ia,j73) (g) is dehned by 

(ii/3i,j2/32|5x (ia '^ ) (p)Ki«i^2a 2 ) 

The basis vectors of Tr in eq. (j6ip are given by 
|ia, i'a'), which belong to the direct product (C x S) x 
(C x S) with dimension 2 4 . It is convenient to arrange 
them as |(cs)J) (v =!,-•• , 16) with 



(62) 



(63) 



(cs)?> = 
Ml> = 

(CB)§) = 
(CB)?) = 

Ml) = 



It, lt> 

U, 1 1) 

2 t, 2 t> 
2 4, 2 t> 
lt,2t) 
(cs)? 1 ) = |U,2t>, 
(cs)? 3 ) = |2t,lt>, 
(cs)? 5 > = |2 4,lt>, 



(cs)|) = |lt,U), 



(cs)|) =|U,U>, 
(c 8 )|> = |2t,2^>, 
(cs)|) = |2|,2|), 
M? > = |lt,2|>, 

M? 2 > = |u,2^>, 
M? 4 > = |2t,U>, 

(cs)? 6 ) = |2 4.,U). 



(64) 



Then, x^ (p, Q) in eq. (|60p can be represented by a 16 x 16 
matrix as 



,(°) = 



- GG FF_ GF — FG 

FF - GG — FG GF 

- GF FG_ GG_ - FF 

FG -GF -FF GG 



(65) 



where each element is a 4 x 4 submatrix given in terms 
of the elements of eq. ([5^)1 ; for example, FF denotes 



FF = 



FffFff FffF)/f FfiFff Ff^Fif 
F n F n F n F u F n F n F n F u 
F it F n FtfFtf F u F n F u F it 



(66) 



with F a/3 F 7 s = F a/3 (p)F lS (p 
r (0) (P,P',9) in eq. (EU) as 

p r (0a) p(06) 

r (o) = 



We can also express 





p(06) p(0a) g 

o o r (0c) 
ooo 







p(0c) 



(67) 



where r (0a - ob ' 0c) are given by 



rW(p,p', g ) = ^^7 6 -^7°, 



(68a) 
(68b) 



p(0c) ( 



p,p ,q) 



i 



'p+p'-q h 

I 5 



2 - 2 
with 1_ denoting the 4x4 unit matrix and 



"1 








r 




"1 








0" 






















1 






















1 








.1 








1_ 













1. 



(69) 



Equation (|63l) may also be transformed into a matrix 
representation. 

Equations ([5(11) and (|61l) form a closed set of self- 
consistent equations for G(p) in the FLEX-S approxi- 
mation, which is applicable to an arbitrary pairing sym- 
metry including a mixture of singlet and triplet pairings. 
By using the representation of eq. (l64l) . practical calcu- 
lations may be performed based on the standard matrix 
algebra plus treatment of integral equations. Some uni- 
tary transformations will be helpful for simplifying those 
calculations. For example, consider the orthogonal ma- 
trix 

~R a -R b 
R b R a 
if 
R c 



R 



(70) 



with 



R° 





r 1 

2 








1-1 

2 







1 
















1 







1 

L 2 








1 




- 1 

2 








1-1 

2 










1 

73 










1 

73 










1 

L 2 








1 

2- 1 



(71a) 



(71b) 
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R c 



1 0" 

75 -75 

7f 71 

1 



(71c) 



By using R, eq. (|57]) is transformed into the diagonal 
form 



f (0) = RT^R- 1 



p(0a) 









ol 


p(06) 


lo 


Ol 


lo 


p(0c) 


Ol 


lo 


lo 


p(0c) 



(72) 



~(0a,06,0c). . 

where 1_ [P,P , <Z) are given by 



p(0a) 



p(0c) 



-ya 








o - 


- 


-v c 














-u c 





. 










v h h 








~-v b 














v c 














-v b 














-v b 



(73a) 

(73b) 
(73c) 



with 



V, 



p+p'-q 



V a (p,p\q) = 2V q - 



v b (p, P ', q )^ Vp+p '-\ Vp - pl . 



v; 



p-p< 



V c (p,p\q) 



Vp-\-p'—q -\- Vp—p' 



(74a) 
(74b) 
(74c) 



The other quantities in eq. (|6"T)) may also be expressed in 
some convenient forms after the transformation. 

Finally, we briefly consider the two-particle Green's 
function defined by eq. (|3"Tj) . which obeys eq. (|4"6")) . Let 
us expand the quantities in eq. (l4*2"j) as eq. (|5"6")) . where 
X^ 0e \p,p',q) is obtained as 

= —Gia,jfl(p)Gj>i3> t ii a >(p-)5p>p 

+G ia<S -j> ^>{^)Gs-j^,i' a '{^-)^' (75) 



Equation (|4"6"]) is thereby transformed into 

pi 

xx i0e) (pi,f,q), (76) 
where T is the Fourier transform of eq. (j45|) obeying 

r (p, p", q) = (p, p',q)-J2 E (ir) (£ Pi > <?1 x (0) (ft , 

pi 

xWi,P',Q)- (77) 



Equations (|76|) and ([77]) correspond to eqs. (6) and (7) 
of Leggett for the spin-singlet pairing in his theory of su- 
perfluid Fermi liquids,— respectively. Here, we can cal- 
culate 1C microscopically based on $ for both the singlet 
and triplet pairings. 



FLEX-S approximation for the contact 
interaction and spin-singlet pairing 



Now, we focus on the contact interaction 



V q = U. 



(78) 



In this case, we can simplify eq. (|61[) for the self-energy 
considerably. Indeed, substituting eq. ((75)) into eq. ([721 . 
we observe that only the 1st, 2nd, 3rd, 4th, 10th, and 
14th diagonal elements are finite in the diagonal matrix 

r ^ ; the relevant 6x6 submatrix is given by 



p(0) 



U 



"1 














0" 





-1 




















-1 




















-1 




















1 




















1 



(79) 



Moreover, absence of momentum dependences in F (0) en- 
ables us to perform the p\ summations in eqs. (j6TJ) and 
([62"T) . Using these properties in eq. ([61"]h we realize that 
we only need the 1st, 2nd, 3rd, 4th, 10th, and 14th rows 
and columns of the two matrices 



~{ia,jf3) / _ 



(p) = RSx 



ip)R- 



(80a) 



(80b) 



These statements on eq. (|6T|) for the contact interaction 
hold true for an arbitrary pairing symmetry. 

We further restrict ourselves to the spin-singlet pairing 
with no magnetic polarization for simplicity. The condi- 
tion implies G aj3 (p) = 5 aP Gp, F n (p) = -F|. t (p), and 
F aa (p) = in eq. ([52")) . It may also be possible for this 
single-component superconductivity to choose the phase 
of F^{p) such that its imaginary part originates solely 
from the ie n dependence as F^i(p) = F^(p*). These 
relations are summarized together with eq. (l53l) as 



G a p(p) = Gf} a {-p) = SapGfi, 

F n (p) = ~F it (p) = F n (p) = -F it (p) = F p = F_f, 
F aa (p)=Q. 

(81a) 

Similarly, the elements of eq. (|54|) will obey 



^ap{p) = T,p a {—p) = 5 Qi aEp-, 

A n (p) = - A it (p) = A n (p) = -A lt (p) = A p = A_jr, 
A aa {p) = 0. 

(81b) 
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Substituting eq. ([81]) into eq. (I50|) . we observe that the 
4x4 matrix equation is divided into a pair of 2 x 2 
equations that are equivalent to one another. The one 
composed of the 1st and 4th rows and columns are given 
by 



Gp 













is,. 



Ep 



-A* 



(82) 



with Gff — G-p and T,p = £_p. 



The other finite subma- 
trix from the 2nd and 3rd rows and columns are obtained 
from eq. (|82l) by changing the signs of Fp and Ap simul- 
taneously. 



Next, we set (ia,j/3) = (If, If), (It) 2 -!-) m e 1- 
to express "Ep and in eq. (|82p in terms of Gp and 
Fp. Elementary calculations with eqs. (j6"3"1) - (j6"6")l . (1701) . 
and (fSTj) show that the relevant 6x6 submatrices of 
eqs. (|80ap and (|80b[) , which are denoted by x^\q) and 
5x^ iot ^ ip)i respectively, are given by 



x (0) 



r (o) 

X- 













-vixf; 





(0) 

xy 




















(0) 

x+ 




















(0) 

x+ 








^xf F 











_y(°) 

^GG 


(0) 

Xff 













(0) 

-Xff 


(0) 

— V - 



(83a) 



UIMT) _ 



&X 



r>(lt,24) 





G 








-G 


V2F 


o~ 







2G 











1 




















4 


-G 








G 


-V2F 







V2F 








-V2F 


-2G 

























F 








-F 


-V2G 


0" 


- 


-2F 





























F 








—F 


-V2G 























_V2G 








-V2G 


2F 


0. 



(83b) 



(83c) 



where the elements of eq. (I83a[) arc defined by 

X [ G ] G (^-Y. G PGP- =X < §U-<1), 
p 

x %W=-J2F f F p _= X ( §U-q), 
p 

p 



XoU-q), 



Xq f (-Q), 



xi%>x^±x^1- 



(84a) 
(84b) 
(84c) 
(84d) 
(84e) 



Let us substitute eqs. (|79|) and (|83|) into eq. (jBTj) for 

{ia,jf3) = (It, It), (It, 2 1). In both cases, the resultant 
Tr can be calculated separately for the two submatrices 
composed of the 2nd, 3rd, and 4th rows and columns 
and the 1st, 5th, and 6th rows and columns. Also using 

J2ffXG%(9) G f-? = Y.qXGG^Gp-qi etc -' in the second 
order, we obtain 

Sp = \un - 2U 2 xf 
q 

+U 2 T,(l^G^+ \^G f X (85a) 



&f= u Y J F.p> - 2U 2 X? ®Fp-t 
p' q 

Here n is the particle density 

n = 2^Gp-e le "°+ 
p 



(85b) 



(86) 



with 0-|- an infinitesimal positive constant^ the func- 
tions Xjf C \ X^\ an( i Xj are defined by 



y(0c) 

—q 



V2 X ( gUq) -x ( gUq) -x { fUq} 



Xq 



0) 



l-U X f{q)' 



xf^xf^i+u^y 1 , 



(87a) 
(87b) 

(87c) 



-q —q v_ —9 

with superscripts c and s denoting "charge" and "spin,'' 
respectively, and Q_ - and Tp are given by 



Gp- V2F ff 6 
V2F ff -2G p 




(87d) 
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5r = 







-V2G P 0' 








(87e) 



respectively. 

Setting F = in eq. (|85al) yields the normal-state self- 
energy 



Si n) - lun 



U 2 



q 



9v (o) | 3 Xgg(^) g p-q 



1 XGG(g) G P-g 

2 1 + ^(9) 



XG G (g)G p -- g - 



i-t/xS® 



(88) 



The second term on the right-hand side completely agrees 
with eq. (2.31) by Bickers and Scalapino in the normal 
FLEX approximation^ as it should. To find the equa- 
tion for T c , let us assume the second-order transition and 
linearize (i) the right-hand side of eq. (|85b[) in terms of 
F and (ii) F itself in eq. (|8"2"j) in terms of A as 



F~ -> F (1) = 
r v ^ r p — 



A p - 



(ie„ 



s£ n) ) 



(89) 



where spin fluctuations are dominant, for example. On 
the other hand, neglecting the particle-hole contribution 
in eq. ([8"5]l altogether by x± (?) ~* and 



.(Oc) 







,«>)_ 
'-GG 







<> -x^Oz) -xg G (£>. 



we obtain the particle-particle ladder approximation 
adopted by Haussmann el al. for the BCS-BEC crossover 
problem^ Hence, eq. ([83)1 enables us to study the 
two cases of considerable physical interest on an equal 
footing, including a crossover from the one limit to 
the other. It is also clear that the FLEX+T-matrix 
approximation; 49 ' 50 which incorporates a contribution of 
the two-particle Green's function into the single-particle 
self-energy, cannot be reproduced as a limit of the FLEX- 
S approximation, where both of the single-particle and 
two-particle Green's functions should be studied start- 
ing from the single functional of eq. (|28p . Hence, a mi- 
croscopic foundation of the FLEX+T-matrix approxima- 
tion remains to be examined critically by including more 
terms from eq. (| 1 T[) beyond the FLEX-S approximation. 



We thereby obtain the linearized self-consistency equa- 
tion as 



IV. SUMMARY 



(i) 



(90) 



where V2g, is given by 



Vfe,=U + U' 



-*>&(?-?)+ 1 XGGif - f) 



2 ^-U X { gUp-p') 



1 X { gUp-p') 
2 ^ + U X { gUp-P') 



(91) 



This V5~, may be regarded as the pair potential for the 
spin-singlet pairing without feedback efiects4^ Note in 
this context that the feedback effects are included natu- 
rally in Eq. (|85b|) . 

Equations (18"2"j) , ([H5]l , and ([83]) form a closed set of self- 
consistent equations for the spin-singlet pairing within 
the FLEX-S approximation. The replacement 



Y (0c) 



X l -\q) 








in eq. (jB"5j) yields the standard FLEX approxima- 
tion for superconductivity with only the particle-hole 
contribution^ -38 ! 40 ' 41 the charge contribution is disre- 
garded further in many cases. This approximation may 
be appropriate for the Hubbard model near half filling 



We have developed a concise self-consistent perturba- 
tion expansion for superconductivity in terms of Green's 
function G and fully symmetrized bare vertex de- 
fined by eqs. © and (|2U| . respectively, so as to incor- 
porate all the pair processes naturally and conveniently. 
Indeed, distinct diagrams for the key functional $[G] of 
eq. (fT7|) are exhausted within the fourth order by the 
five graphs of Fig. [5] Corresponding analytic expressions 
can be written down immediately as eq. (|26l) based on 
the Feynman rules (a)-(f) given around eq. (|2"3"|) . Using 
an approximate $ thereby constructed, we can calculate 
G self-consistently by eqs. (O and (TT3")) . the thermody- 
namic potential subsequently by eq. (jlip . and also the 
two-particle Green's function of eq. (|3Tj) by eq. (|44|) with 
eqs. (|25b|) . (f3T)|). and flU]) . Moreover, entropy may be 
written down explicitly once G and £ are knowni^ 9 . Thus, 
a single procedure of choosing an approximate $ enables 
us to calculate the whole thermodynamic hierarchy, and 
the results may be improved by starting from a better 
$. Various fluctuations and competing orders may be 
described simultaneously through the matrix structure 
of r*- ). Since the conservation laws are obeyed auto- 
matically, the basic formalism of §2 in the coordinate 
representation can be applied to study non-equilibrium 
phenomena by changing the Matsubara contour into the 
Schwingcr-Kcldysh contour* 1 ^. The idea of symmetrizing 
the bare vertex in terms of the particle-hole indices will 
also be useful for the electron-phonon Hamiltonian to 
simplify its perturbation expansion. 

This formalism has been transformed for homoge- 
neous systems into the momentum-energy representa- 
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tion in §3.1. We have subsequently derived a closed 
set of self-consistent equations for the contact poten- 
tial in the spin-singlet FLEX-S approximation as eqs. 
([82]). (155]) . and which adds all the pair processes 

characteristic of superconductivity to the normal FLEX 
approximation.— As shown in the paragraph after eq. 
(|9ip. it contains extra terms besides those in the stan- 
dard FLEX approximation for superconductivity, and 
also embraces the particle-particle ladder approximation 
by Haussmann et al^ for the BCS-BEC crossover prob- 
lem. Thus, our FLEX-S approximation incorporates the 
particle-hole and particle-particle diagrams on an equal 
footing with their mixing effects, thereby enabling us to 
study the BCS-BEC problem from low- to high-density 
regions continuously. 
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Appendix A: Perturbation expansion with 

The Feynman rules for the normal-state expansion of 
eq. (|17l) with eq. (|15l) are summarized as follows. 

(a) Draw all the nth-order closed skeleton diagrams 
that are topologically distinct. For each such di- 
agram, associate the factor (— l) ,l+1 /4 n n!/3. 

(b) For each vertex, associate r(°)(AA', vv'). 

(c) Identify the number C na of possible connections of 
vertices for diagram a under consideration. 

(d) Consider a specific connection of vertices for dia- 
gram a and associate G(A', fj) for each line arriving 
at A' from fj. 

(e) Identify the number l na of permutations to real- 
ize the connection. We may find it by replac- 
ing r(°)(lT,22') ->■ 6(\,\')6(2,2')V(r x - r 2 ) and 
equating l na with the number of closed particle 
loops in the resultant diagram with V ^ 

(f) Multiply the expression by (-l)'™C, lo . 

These rules can also be used for superconductors, where 
the number £ na for an "anomalous" diagram may only 
be found by a direct counting of permutations. 

A key quantity in the rules is the numerical weight fac- 
tor C na - Since the two operators ifiiipi in eq. (fT"5|) with 
the same subscript i (— 1,2) are equivalent, the factor 
may be obtained elementarily by a combinatorial consid- 
eration. For example, the factors for Fig.[3][a) and (b) are 
trivially identified as C\ a = 2 and C\ b — 1, respectively. 
Together with l\ a = 2 and l\ b — in terms of eq. (TTT)1) . 
we arrive at eq. (TT5)) for Fig. [3] following the above rules. 



Similarly, the weight factors for Fig. [T] (a)-(j) are found 
as C 3a = 2 5 2!, C 3b = 2 3 2!, C 3c = 3(2 3 ) 2 , C 3d = 3!(2 3 ) 2 , 
C 3e = 3!25, C 3f = 3(2 3 ) 2 , C 3g = 3!(2 3 ) 2 , C 3h = 3!2 3 , 
C 3 i — 3 • 2 5 , and C 3 j — 2 6 , each of which corresponds to 
the number of possible connections of three vertices for 
realizing the relevant topological structure. We thereby 
obtain the relative weights of Fig. [TJ 

Appendix B: Momentum-energy representation for 
lattice models 

Here, we discuss the momentum- energy representation 
of homogeneous lattice models. First of all, momentum 
integrations in eqs. (H5J and (|55l) should be restricted to 
the first Brillouin zone. This fact may be incorporated 
manifestly in these formulas by introducing the function 

e p ^/ J^s 3 (k- P ) (bi) 

on their right-hand sides. 

With these trivial modifications, we consider the basic 
two-particle quantities in eq. (|25[) . Substituting eq. (|16[) 
into eq. (|20b|) . we can write the matrix element of eq. 
(|25a[) explicitly in terms of V as 

(ii^|r(°)|22; y ) 

= ~[(-l) <+ %^(ri -r' 2 )S(r[ - r' 2 )6(l,l')6(2' ,2) 

-(-iy+%-^(r 2 - r' 2 )S(r 2 - 4)6(1, 2)6(2', 1') 
+(-l) i+i '6 j ., 3 - i 6 j , 3 - i .V(r' 2 ~ r[)6(^ - t{)6(1,2') 
x6(l',2)]. (B2) 

Let us substitute eqs. (|48cl) and (|55|) with eq. (|B1[) into 
eq. (|B2[) and transform the resultant expression into the 
form of eq. (|56p by changing integration variables. We 
thereby obtain the Fourier coefficient in eq. (I56|) as 

(ia,i'a'\TW(p,p',q)\jP,j'l3') 

— (— l) l+l 6ij6i> j> 5 a p6 a > p'QpOp-qQp-p'Vp-p' 

+(— l) l+t 6j> , 3 -i6j y3 -i'S a p'6 a 'i3QpQp'Qp + p'^q 

xW-J- ( B3 ) 

Note that p 1 and q above may lie outside the first Bril- 
louin zone to give a finite value for the matrix element. 
Similarly, the Fourier coefficient of eq. (|59[) for lattice 
models is obtained as 

(ia,i'a'\x (0) (p,q)\jPJ'P') = -G ia ^(p)G r ^, ila ,(p^) 

xQpOp-q. (B4) 

By incorporating functions accordingly in other for- 
mulas, the whole results of §3 also become applicable to 
lattice models. 
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